The nuclear quadrupole moment of 57 Fe from microscopic nuclear and atomic 

calculations 
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The nuclear quadrupole moment of the I n —3/2~ excited nuclear state of 57 Fe at 14.41 keV, impor- 
tant in Mossbauer spectroscopy, is determined from the large-scale nuclear shell-model calculations 
for 57 Fe and also from the electronic ab initio and density functional theory calculations including 
solid state and electron correlation effects for the molecules Fe(CO)s and Fe(C5Hs)2. Both inde- 
pendent methods yield very similar results. The recommended value is 0.16(1) eb. The NQM of the 
isomeric 10 + in 54 Fe has also been calculated. The new value (0.5 eb), consistent with the perturbed 
angular distribution data, is by a factor of two larger than the currently recommended value. 
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Mossbauer spectroscopy of 57 Fe plays an important 
role in the structural determination of iron contain- 
ing solid state compounds. In principle, the nuclear 
quadrupole moment (NQM) of the isomeric 7=3/2 state 
in Fe can be determined from Mossbauer data; how- 
ever, the analysis requires the calculation of the elec- 
tric field gradient (EFG). As these atomic calculations 
are quite involved, studies employing different methods 
arrived at quite distinct results; the values of NQM in 
the range from -0.19 to +0.44 eb have been reported |Q. 
This quite unsatisfying situation could also not have been 
settled by nuclear structure calculations of the NQM as 
calculations within the nuclear shell model, the most reli- 
able tool for such studies, had to be performed in strongly 
truncated model spaces and with rather untested effec- 
tive interactions. In recent years, decisive progress has 
been achieved in both the atomic calculations of the EFG 
and in nuclear shell model studies. 

In 1995, Dufek and co-workers applied the density 
functional theory for a series of iron-containing solid state 
compounds. For 57 Fe they obtained a NQM of 0.16 eb 
|^| , in contradiction with the previously accepted value of 
0.082 eb obtained from Hartree-Fock (HF) EFG calcula- 
tions H and from truncated nuclear shell-model calcula- 
tions combined with the perturbed angular distribution 
data ||. In subsequent work, Su and Coppens obtained 
a NQM value of 0.12(3) eb using Sternheimer-corrected 
EFGs ||. In this Letter, we shall demonstrate that state- 
of-the-art nuclear and atomic physics calculations lead 
to the same NQM for the 57 Fe isomeric state, settling a 
long-standing controversial issue. 

We shall begin with the nuclear physics discussion. 
The nucleus 57 Fe has two low- lying 3/2~ states which 
are experimentally split by only 353 keV. To describe 



the structure of such nearly degenerate states is quite 
demanding. Clearly, the large-scale shell model is the 
method of choice. Due to recent progress in programming 
and hardware development, modern shell- model calcu- 
lations based on microscopic effective interactions can 
handle configuration spaces that were prohibitively large 
only several years ago ^]. More specifically, modern di- 
agonalization shell-model codes can now handle medium- 
mass nuclei (^4=50-60) in the middle of the pf shell in 
full OTiuj space. To put things in perspective, the shell- 
model calculations of the Utrecht Group M , which were 
used in Ref. Q to extract the NQM of 57 Fe, restricted 
the number of holes in the fa/ 2 orbit to three. 

Shell-model calculations depend crucially on two fac- 
tors: the model space and the effective interaction. 
Our calculations have been performed using the code 
NATHAN g. NATHAN has been developed in the jj- 
coupling scheme using quasispin formalism. We adopt 
a version of the code adapted to shared-memory parallel 
machines. For the two lowest 3/2~ levels, we assumed the 
truncated space (containing 8,120,105 J 7r =3/2~ states) 
in which maximally 6 nucleons were allowed to be ex- 
cited from the fa / 2 orbital to the rest of the pf shell. To 
test the convergence of our results, we have performed a 
full pf shell calculation for the lowest 3/2~ state. This 
calculation includes 25,743,302 states (^2-10 13 non-zero 
matrix elements) and is one of the largest shell-model 
diagonalization performed to date. For all adopted inter- 
actions, the results of the truncated and complete calcu- 
lations are identical. 

During the last few years, a considerable effort went 
into the development of effective interactions in the pf 
shell. In this work three different effective interactions 
have been employed: KB3F, KB3G, and FDP6. The in- 
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teractions KB3F || and KB3G [|Tc]] are both reasonable 
attempts to correct the defects of the well-known KB3 
interaction |Tl| ] in the upper part of the p/-shell. They 
differ slightly in their collectivity of states around 58 Ni. 
As far as 57 Fe is concerned, KB3G is clearly the interac- 
tion of choice; it has been shown to be very successful in 
describing experimental data (including energy levels and 
electromagnetic properties) in the mass region A=50-52 
and around A=56 0. The FDP6 interaction jfj] was 
originally fitted to the spectroscopic properties of 
nuclei and has since been extended to nuclei around the 
N=Z=28 shell closure |13|. For the effective charges, 



we took the quadrupole charges 1.5e for protons and 
0.5e for neutrons when calculating E2 transitions and 
moments, and the spin and orbital gyromagnetic factors 
5 s =0.755g arc , gl=l.l Mat, and g l u =-0.1 ^ N for Ml tran- 
sitions and moments, e.g. E3|. We adopt the oscillator 
parameter b = 1.01 A 1 / 6 fm. 

Our shell-model calculations reproduce well the exper- 
imental level scheme of 57 Fe. In particular, the low-lying 
7^=3/2^ doublet is predicted by theory. Experimentally, 
the energy difference between these states is 353 keV, 
while theoretically it is 526, 209, and 465 keV using the 
KB3F, KB3G, and FPD6 interactions, respectively. The 
third 3/2 _ state appears much higher in energy, at about 
1.6 MeV. Table | summarizes our shell- model results for 
the quadrupole and magnetic moments of the two lowest 
3/2 ~~ states in 57 Fe. It is seen that these close- lying states 
in the doublet are predicted to have very different shapes 
and magnetic moments, i.e., their E2 and Ml moments 
have different signs. Consequently, even small changes 
in the shell-model interaction, hence in the coupling be- 
tween these states, significantly impact theoretical pre- 
dictions. We note that a convergence of the shell-model 
results is not reached if only three f?/2 holes are allowed, 
as was done in the previous truncated calculations 0. 
For example, the value of the NQM of the 3/2 j~ state, 
Q(3/27), obtained in such truncated calculations with 
the KB3G interaction, is 0.097 eb, as compared to the 
exact value of 0.064 eb. 

In order to assess the sensitivity of the results on the 
adopted effective force, we computed overlaps between 
wave functions of the 3/2~ states obtained with the dif- 
ferent interactions. It turned out that the spaces of the 
two lowest 3/2~ states in KB3F and KB3G are practi- 
cally (up to 95%) the same. That is, both 3/2~ wave 
functions obtained with KB3F can be approximately de- 
rived from the 3/2~ doublet calculated with KB3G by 
simple rotation. This does not hold for KB3G and FDP6. 
Since FPD6 tends to excite particles to the f 5 / 2 shell 
rather than the p^/ 2 shell (FPD6 puts the f 5 / 2 orbit too 
low at N—28), the contributions from higher-lying 3/2~ 
states amount to about 37%. 

Due to the great sensitivity of the shell-model predic- 
tions caused by the not-very- well-controlled off-diagonal 
coupling between the 3/2 _ doublet, the calculations need 



to be constrained by the available experimental data. To 
this end, we choose the magnetic moment of the 3/2^ 
state, which has been precisely determined experimen- 
tally, ^(3/27)^-0.1549(2) fj, N Assuming that the "true" 
wave functions of the 3/2 _ doublet are given by a simple 
rotation of the shell-model states, 
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one can determine the mixing amplitudes a% and a 2 
{a\ + a 2 = 1) by requesting that the measured value 
of //(3/2T/) be reproduced. The results of such two-level 
mixing calculations are displayed in Table |j[ Contrary 
to the pure shell-model results, one obtains a satisfying 
agreement between KB3F and KB3G, and the data. In 
particular, the quadrupole moments predicted by these 
interactions are very similar, (5(3/2^)^0.17 eb. The only 
serious difference is the sign of the magnetic moment of 
the 3/27 state. Unfortunately, this quantity is poorly 
determined experimentally: /i(3/27)<0.6 jtxjv The predic- 
tions of FPD6 for the E2 moments are rather far from 
the data: both 3/2 _ states are calculated to be prac- 
tically spherical, and the B(E2) transition connecting 
these states is enhanced by a factor of ~40. Part of this 
failure comes from too-low a placement of the / 5 / 2 level 
in FPD6, as discussed in detail in Ref. |Tc|] , where more 
examples of KB3G and FPD6 calculations can be found. 
Therefore, we conclude that the best shell-model esti- 
mate of the quadrupole moment of the 3/27/ state in 57 Fe 
is 0.16 eb. 

In Ref. jl6| it has been pointed out that the ratio of 
quadrupole moments of the isomeric 10+ state in 54 Fe to 
the 3/2- state in 57 Fe, Q[ 54 Fe(10+)]/Q[ 57 Fe(3/27)], rep- 
resents a stringent constraint. Possible uncertainties in 
the atomic EFG calculations cancel largely in the exper- 
imental determination of the ratio from Mossbauer data. 
The most recent value is Q[ 54 Fe(10+)]/Q[ 57 Fe(3/27)] 
=3.62± 0.22 0. Unfortunately, the NQM of the iso- 
meric state in Fe is not known experimentally. We 
have thus calculated the NQM of the isomeric 10 + state 
in 54 Fe in a large-scale shell-model study in which we al- 
lowed 8 particles to be promoted from the fj/ 2 orbital 
to the rest of the shell. We have checked the conver- 
gence of our results by performing calculations in which 
only 6 particles could be promoted; the NQMs were the 
same at both levels of truncation. For all the interac- 
tions used, we find rather similar NQMs: 0.51 eb (KB3F), 
0.50 eb (KB3G), and 0.56 eb (FDP6). Using the NQMs 
of Table |l| for the isomeric state in 57 Fe, these values 
correspond to Q[ 54 Fe(10+)]/Q[ 57 Fe(3/27)] ratios which 
agree rather well with the measurement for the KB3F 
and KB3G interactions (2.9 and 3.1, respectively). For 
the FDP6 interaction, the ratio is significantly too large, 
as the NQM of the isomer in 57 Fe is severely underesti- 
mated. We note that our predicted values for the mag- 
netic moment of the 10 + state in 54 Fe (6.4 /z^r in KB3F, 
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6.5 hm in KB3G, and 7.1 /ijv in FDP6) compare well with 
the experimental value of 7.28±0.01 n.m fill. 

We now turn to our atomic physics studies. The use 
of the density functional theory for the calculation of 
EFGs in transition metals is questionable. Recent calcu- 
lations on CuF showed a variation of the Cu EFG rang- 
ing from +0.50 au for the local density approximation 
(LDA), and +0.44 au at the generalized gradient approx- 
imation (GGA) level (BPW91) to +0.07 au at the hybrid 
level (B3PW91) |lj],|2(]], as compared to the experimen- 
tal value of —0.31(2) au pi"| . In contrast, relativistic ab 
initio coupled-cluster calculations give —0.34 au |19|, in 
perfect agreement with the experimental result. 

In this work, we carried out density functional as well 
as ab initio calculations for the molecules Fe(CO)s and 
Fe(C5H 5 ) 2 . We have adopted a wide range of exchange 
and correlation functionals for the electronic structure 
calculations of the free molecules Fe(CO)5 and Fe(C5H 5 ) 2 
(for the terminology see Refs. [|l9|,^2)): Hartree-Fock- 
Slater (HFS), Xa, Local Density Approximation (LDA), 
the GGA functionals B-HFS (Slater exchange plus Becke 
nonlocal exchange), B-LYP (B-HFS plus the correlation 
functional of Lee, Yang and Parr), B-PW91 (B-HFS plus 
the correlation functional of Perdew and Wang), the hy- 
brid GGA functionals B3-LYP (Becke three-parameter 
functional), B3-PW91 (same as B3-LYP except with the 
nonlocal correlation term of Perdew and Wang), BHH 
(Becke half-and-half together with the LYP correlation 
functional), and BHH- LYP (same as BHH but with 0.5 
of the Becke nonlocal exchange term added to the en- 
ergy). For comparison with DFT, we carried out the HF 
many-body perturbation theory (MBPTn) up to order 
n=4, as well as coupled-cluster singles doubles including 
non-iterative triple excitations [CCSD(T)] pSj ]. The elec- 
tronic coupled cluster calculations required three months 
of CPU time and 20 Gbytes of disk storage on a Origin 
2000 SGI. For Fe(CO) 5 we investigated solid state ef- 
fects to the iron electric field gradient by performing HF, 
B3LYP, and LDA calculations using the program CRYS- 
TAL98 [^4| and the solid state structure given by Bose 
and Blaser (29]. Detailed structural data and basis sets 
used will be published elsewhere ^(J . 

The results of our electronic structure calculations for 
the iron EFGs of Fe(CO)5 and Fe(C5H 5 ) 2 are shown 
in Table 



[II 



We first note that the single-reference 
many-body perturbation theory shows extreme oscilla- 
tory behavior and is practically useless for the determi- 
nation of transition element EFGs. DFT results range 
from 1.09 au to 1.57 au for Fe(CO)5, and from 1.36 au to 
2.49 au for Fe(C5Hs)2, depending on the functional ap- 
plied. If we accept the coupled-cluster result as the most 
accurate value, the hybrid GGA functionals perform well 
for Fe(CO)5, while the non-hybrid GGA functionals are 
preferred for Fe(C 5 H 5 ) 2 . We note that for CuCl the BHH 
functionals gave the best description |H| . This is clearly 
not a satisfying situation. On the other hand, an en- 



couraging result is that solid state effects from nearest- 
neighbor interactions in Fe(CO)s are small and can basi- 
cally be neglected. Based on the CCSD(T) EFGs, we ob- 
tain from the Mossbauer data of Fe(CO)5) (2.57 mm/sec 
|3) andFe(C 5 H 5 ) 2 (2.4 mm/sec Q) aNQM of 0.177 eb 
and 0.159 eb, respectively. However, basis set incomplete- 
ness and relativistic effects may increase the iron EFGs 
and, therefore, further decrease the NQM | p6[ . Con- 
sequently, our best estimate using EFGs together with 
Mossbauer data for these molecules is 0.15(2) eb. In or- 
der to improve further on these results, large-scale rela- 
tivistic coupled-cluster calculations are necessary, which 
are currently not feasible for such big molecules. We em- 
phasize, however, that EFGs obtained from current DFT 
for transition element compounds should be taken with 
some care as the results in Table III show. 



In summary, the quadrupole moment of the first 3/2 _ 
state in 57 Fe at 14.41 keV, important in Mossbauer spec- 
troscopy, has been calculated using state-of-the-art nu- 
clear and atomic models. Both calculations yield re- 
sults which are consistent with Q(3/2j~)«0.16eb, in nice 
agreement with the recently reported value [^) and the re- 
sults of Ref . H . As a by-product of this work, the NQM 
of the isomeric 10 + in 54 Fe has also been calculated. The 
new value, consistent with the perturbed angular distri- 
bution data, is 0.5 eb, i.e., a factor of two greater than 
the currently recommended value |17|. 

This research was supported by the U.S. Department 
of Energy under Contract Nos. DE-FG02-96ER40963 
(University of Tennessee), DE-FG02-97ER41019 (Uni- 
versity of North Carolina), DE-AC05-00OR22725 with 
UT-Battcllc, LLC (Oak Ridge National Laboratory), 
DE-FG02-92ER40694 (Tennessee Technological Univer- 
sity), the Danish Research Foundation, the Schweizer 
Nationalfond (grant 20-61822.00), the Marsden Fund 
Wellington, the Royal Society of New Zealand, and the 
Auckland University Research Committee. The Auck- 
land group is indebted to the Alexander von Humboldt 
Foundation for financial support. GMP wants to thank 
the Carlsberg Foundation for a fellowship. 



3 



V.S. Rusakov and D.A. Khramov, Bull. Russ. Acad. Sci. 
Phys. 56, 1118 (1992). 

P. Dufek, P. Blaha, and K. Schwarz, Phys. Rev. Lett. 75, 
3545 (1995). 

K.J. Duff, K.C. Mishra, and T.P. Das, Phys. Rev. Lett. 
46, 1611 (1981). 

S. Vajda, G.D. Sprouse, M.H. Rafailovich and J.W. Noe, 

Phys. Rev. Lett. 47, 1230 (1981). 

Z. Su and P. Coppens, Acta Cryst. A 52, 748 (1996). 

K. Langanke and A. Poves, Nucl. Phys. News 10, 3 

(2000). 

R. Vennink and P.W.M. Glaudemans, Z. Phys. A 294, 
241 (1980). 

E. Caurier, G. Martmez-Pinedo, F. Nowacki, A. Poves, 
J. Retamosa, and A. P. Zuker, Phys. Rev. C 59, 2033 
(1999). 

H. Kaiser et ai, Nucl. Phys. A669, 368 (2000). 

A. Poves, J. Sanc hez-Solano, E. Caurier, and F. Nowacki, 



nucl-th/0012077| ; to be published. 

G. Martmez-Pinedo, A. P. Zuker, A. Poves, and E. Cau- 
rier, Phys. Rev. C 55, 187 (1997). 

W.A. Richter, M.G. Vandermerve, R.E. Julies, and B.A. 
Brown, Nucl. Phys. A523, 325 (1991). 
T. Otsuka, M. Honma, and T. Mizusaki, Phys. Rev. Lett. 
82, 1588 (1998). 

A. Bohr and B.R. Mottelson, Nuclear Structure, Vol. II 

(Benjamin, New York, 1975). 

M.R. Bhat, Nucl. Data Sheets 85, 415 (1998). 

S. Vajda, G.D. Sprouse, M.H. Rafailovich, and J.W. Noe, 

Phys. Rev. Lett. 47, 1230 (1981). 

M. Hass et ai, Nucl. Phys. A414, 316 (1984). 

M.H. Rafailovich, E. Dafni, J.M. Brennan and G.D. 

Sprouse, Phys. Rev. C 27, 602 (1983). 

P.Schwerdtfeger, J. K. Laerdahl, and M. Pernpointner, 

J. Chem. Phys 111, 3357 (1999). 

E. van Lenthe and E. J. Baerends, J. Chem. Phys 112, 
8279 (2000). 

J. Hoeft, F.J. Lovas, E. Tiemann, and T. T0rring, Z. 
Naturforsch 26a, 240 (1970). 

Recent Developments and Applications of Modern Den- 
sity Functional Theory, Theoretical and Computational 
Chemistry, ed. by J.M. Seminario (Elsevier, Amsterdam, 
1996). 

M.J. Frisch et al., program Gaussian 98, Gaussian Inc., 
Pittsburgh PA, 1998. 

V.R. Saunders et al., program CRYSTAL98, University 
of Torino, Torino (1998). 

R. Bose and D. Blaser, Z. Krist. 193, 289 (1990). 

P. Schwerdtfeger, T. Sohnel, M. Pernpointner, and J. K. 

Laerdahl, to be published. 

R. Greatrex and N.N. Greenwood, Discuss. Faraday Soc. 
47, 126 (1969). 

R.L. Collins, J. Chem. Phys. 42, 1072 (1964). 
M. Fanciulli et al, Phys. Rev. B 59, 3675 (1999). 



TABLE I. Results of the shell-model calculations for the 
quadrupole and magnetic moments of the two lowest 3/2~ 
states in 57 Fe using effective interactions KB3F, KB3G, and 
FPD6. The experimental value of //(3/2T/) is -0.1549(2) /xjv 
fl 







KB3F 


KB3G 


FPD6 
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3/2r 


0.16 


0.06 


-0.17 


(eb) 


3/2 2 - 


-0.16 


-0.07 


0.17 




3/2- 


-0.32 


-0.49 


-0.51 


(MJv) 


3/2" 


0.26 


0.23 


0.49 



TABLE II. Results of the two-level mixing calculations 
(mixing amplitudes, quadrupole and magnetic moments, and 
transition probabilities) for the two lowest 3/2 _ states in 
57 Fe using effective interactions KB3F, KB3G, and FPD6. 
The mixing amplitudes ot\ and a?2 (0) for the 3/27/ state 
have been adjusted to reproduce the experimental value of 
p(3/2i )=-0.155 hn B{M1) and B(E2) denote 3/27 -> 3/27/ 
transition rates. See text for details. 
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0.07 5 



TABLE III. Calculated iron electric field gradient for 
Fe(CO)s and Fe(CsHs)2 at various levels of theory. SS de- 
notes solid state calculations. All values are in au. 
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Fe(C 5 H 5 )2 
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HF 


1.394 


3.157 


SSHF 


1.367 




MBPT2 
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MBPT3 


-0.740 


3.385 


MBPT4 


9.484 


-0.250 


CCSD(T) 


1.394 


1.496 






DFT 


Xq 


1.148 


1.374 


HFS 
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1.363 


LDA 


1.154 
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1.122 




BHFS 


1.187 


1.434 


BLYP 


1.203 


1.425 


BPW91 


1.203 


1.377 


B3LYP 


1.393 


1.854 


SSB3LYP 


1.357 




B3PW91 


1.395 


1.806 


BHH 


1.547 


2.429 


BHHLYP 


1.573 


2.488 
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